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In contrast to perturbative QCD, the analytic QCD models have running coupling whose analytic 
properties correctly mirror those of spacelike observables. The discontinuity (spectral) function of 
such running coupling is expected to agree with the perturbative case at large timelike momenta; 
however, at low timelike momenta it is not known. In the latter regime, we parametrize the unknown 
behavior of the spectral function as a sum of (two) delta functions; while the onset of the perturbative 
behavior of the spectral function is set to be 1.0-1.5 GeV. This is in close analogy with the "minimal 
hadronic ansatz" used in the literature for modeling spectral functions of correlators. For the running 
coupling itself, we impose the condition that it basically merges with the perturbative coupling at 
high spacelike momenta. In addition, we require that the well-measured nonstrange semihadronic 
(V + A) tau decay ratio value be reproduced by the model. We thus obtain a QCD framework 
which is basically indistinguishable from perturbative QCD at high momenta (Q > 1 GeV), and at 
low momenta it respects the basic analyticity properties of spacelike observables as dictated by the 
general principles of the local quantum field theories. 



o 



X 



PACS numbers: 12.38.Cy, 12.38.Aw,12.40.Vv 



I. INTRODUCTION 
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f^ ' The general principles of the local quantum field theories imply [J 0] that the spacelike observables V(Q 2 ) are 

analytic functions of Q 2 in the complex Q 2 -plane, with the nonanalyticity allowed only on the timelike semiaxis Q 2 < 0. 
Here, — q 2 = Q 2 where q is the typical momentum of the considered proccess, e.g., momentum of the exchanged photon 
in deep inelastic scattering. However, in perturbative QCD (pQCD), the running coupling a pt (<2 2 ) = a s (Q 2 )/ir in 
p^ ' most renormalization schemes has (Landau) singularities in the complex Q 2 -plane close to the origin (\Q 2 \ < 1 GeV 2 ) 

^- , which do not reflect the aforementioned analytic properties. This fact represents a serious, and often underestimated, 
conceptual and practical problem. The (spacelike) observables evaluated in pQCD, as truncated power series of 
flpt(tQ 2 ) (with k ~ 1), do not possess the aformentioned analyticity properties; furthermore, for low values of 
\Q 2 \ (~ 1 GeV 2 ), they are numerically unreliable due to the vicinity of the scale kQ 2 to the unphysical (Landau) 
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singularities. This numerical unreliability reflects itself in a strong dependence on the renormalization scale (-H- n) 
C*~) ■ and scheme. 

On the other hand, studies using Dyson-Schwinger equations [3, 4] and lattice calculations 5] indicate that the QCD 
running coupling is finite ( "conformal" ) at Q 2 = 0, and has, at least at positive Q 2 values, no Landau singularities. 

The problem of Landau poles in the QCD coupling was first addressed in a systematic manner about 15 years ago 
by the authors of p-Q, who constructed and used an analytic QCD coupling parameter A\ (Q 2 ) closely based 
on the perturbative coupling parameter a pt (Q 2 ): in the dispersive integral expression for a pt (Q 2 ), they removed the 
integration over the offending spacelike discontinuity cut (i.e., at —Q 2 = a < 0), while keeping the discontinuity 

(spectral) function p^ (a) = Im a pt (Q 2 = —a — ie) unchanged on the timelike discontinuity cut, i.e., at a > 0. 
Therefore, we can call this the Minimal Analytic (MA) model. The authors of \Q%] called it the Analytic Perturbation 

Approach (APT), and provided an analogous method of construction of analytic analogs A n (Q 2 ) of the powers 

Opt(Q 2 )™ for n = 2, 3, A method of construction of analytic analogs A v (Q 2 ) of noninteger powers a p t{Q 2 ) u for 

MA (APT) was developed by the authors of Q . Yet another analytic QCD model, based on the minimal analytization 
of the beta function da p t(Q 2 )/d\nQ 2 , was constructed and used in flOf . 

On the other hand, a more general approach of constructing analytic QCD coupling Ai(Q 2 ), based on a given 
choice of the discontinuity function pi(a) = Im Ai(Q 2 = —a — ie) (for a > 0) was emphasized in [IJ QjJ and (in a 
more specific context) in [lj|]. The spacelike coupling A\{Q 2 ) is then constructed from pi(cr) by the usual dispersion 
relation 
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Dedicated to Olivier Espinosa (1961-2010). 



From our standpoint it is this analytic coupling Ai (or equivalently: the spectral function pi) that defines the analytic 
QCD (anQCD) model. The construction of analytic analogs A n (Q 2 ) of integer powers a p t(Q 2 ) n , applicable to any 
such analytic QCD model, was performed in [TUlL^I], usm g the relations between the logarithmic derivatives a p t,n(Q 2 ) 
(ex d"- 1 a pt (Q 2 )/d(lnQ 2 )™^ 1 ) and powers a pt (Q 2 ) k 






a™ t = a pt ,„ + 2_^ k m (n)a vt ,n+m , (2) 

m>l 

and 1 imposing the condition of analogy on the logarithmic derivatives of a pt and of A\. a p t.n+m l— ^ ■A-n+m- This 
condition was shown to be imperative, in order to keep the renormalization scale and scheme dependence of the 
resulting truncated analytic series (in terms of A n , or A n ) of physical observables under control. The construction of 
the power analogs A n of a" t in general analytic QCD models was obtained thus from A\ via the relations analogous 
to © 2 

m>l 

The extension of this construction to noninteger power analogs A U {Q 2 ), for general analytic QCD models, was 
performed in [17| . 

So, from our standpoint, it remains an outstanding problem to obtain or construct the most acceptable analytic 
coupling A\{Q 2 ), or equivalently, the spectral function pi(<r) = Im Ai(Q 2 = — er — ie) (for a > 0). It is reasonable to 
assume that at large a (> 1 GeV 2 ) we have pi(cr) = Pi (c), i.e., the spectral function agrees with the pQCD result. 
On the other hand, at low a < 1 GeV 2 , the exact behavior of p\ (a) is unknown. 

The construction of A\{Q 2 ) can be performed in two different ways. One way is to construct first the beta 
function f3(Ai) — dAi / d In Q 2 as function of A\. This approach is convenient if we take the position that j3(x) is 
an analytic function of x at x = 0. In such a case, it tuns out that we obtain perturbative QCD, i.e., A n = A\. 
However, in this case, after ensuring additionally the analyticity of A\(Q 2 ) as a function of Q 2 , it turns out to be very 
difficult to reproduce the correct measured value of the tau lepton (nonstrange) semihadronic (V + A) decay ratio 
r T sa 0.20, cf. [18]. In fact, as shown in [18| . the large enough value r T can be obtained in the perturbative analytic 
QCD frameworks only at an (unacceptable?) price of choosing a renormalization scheme with increasingly large (3j 
coefficients, which makes the analytic perturbation series of observables convergent only when up to four terms are 
included, and the fifth (N 4 LO) term in the expansion shows an explosive increase. 

Another way is to construct first the discontinuity (spectral) function pi(<r) = lmAi(Q 2 = — a — ie) (for a > 0). 
This approach leads, in general, to nonperturbative analytic QCD, i.e., A n turns out to be different from A\. We can 
follow here analogous ideas used in the construction of the spectral functions of spacelike observables (correlators) in 
the literature, e.g., |19l - [22l |. In these references, analytization is applied directly to a considered (spacelike) observable 
T>(Q 2 ) itself. Some of the new nonperturbative parameters introduced there were thus specific to the chosen observable. 
On the other hand, we take here the standpoint that it is the (universal) QCD coupling that needs analytization; 
and that the additional nonperturbative contributions for a considered observable, not contained in the analytized 
leading-twist contribution, are accounted for by a procedure containing other universal parameters. Such parameters 
can be vacuum expectation values of higher dimensional operators, and the aforementioned additional nonperturbative 
contributions are represented by higher- twist terms of the Operator Product Expansion (OPE) |23| . 

In our approach, at a large enough threshold value Co = Mq (~ 1 GeV ) we have the onset of the perturbative 
behavior for the discontinuity function of the coupling 

Pl (a)=p[ pt \a), ior<r>M 2 . (4) 

On the other hand, in the regime < a < Mq, the behavior is nonperturbative and unknown in detail, and could be 
parametrized as a sum (with different weights) of delta functions 

n 
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pi (a) = > F 2 5(a - M 2 ) , for < a < M 2 . (5) 
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The recurrence relations leading to the above relations, within the context of the MA (APT) model of [34al! were given in 1141, 11511 . 
The construction of higher power analogs An as linear operations on Ai (not as: A±) incorporates a nice functional property: its 
compatibility with linear integral transformations, such as Fourier or Laplace [lq| . 



As has been argued in [2J, |25[ , introduction of a sufficient number of positive delta functions in the discontinuity func- 
tion p(a) = Im/(Q 2 = — a — ie) can approximate sufficiently well any positive Stieltjes function f{Q 2 ). The analytic 
coupling Ai(Q 2 ) is a positive Stieltjes function [231. On the other hand, application of such type of approximations 
has been applied to the spectral functions of certain current correlators |2l| , under the name of the "minimal hadronic 
ansatz." 

In [26j |. we constructed in this way a simple one-delta analytic QCD model, by introducing one delta function in 
the low-cr regime. The model contains three free parameters, 3 which were fixed by the condition (two requirements) 
A\ — a pt ~ (A 2 /Q 2 ) 3 at large \Q 2 \ > A 2 (where A 2 = Aq CD ~ 0.1 GeV 2 ), and the (one) requirement of reproducing 
the correct value of the semihadronic tau decay ratio r T . 

It may be regarded as overly optimistic to approximate the unknown low-cr regime by a single delta funtion. In the 
present work, we go beyond the one-delta approximation, and investigate how to parametrize the low-cr regime in a 
reasonable manner with two positive delta functions. Since such an extension introduces in the model several new 
parameters, we will fix those parameters by specific reasonable conditions, which will be similar in spirit to those of 
the one-delta case. Specifically, the model can be made even closer to pQCD, A\ — a pt ~ (A 2 /Q 2 ) 5 at large \Q 2 \ > A 2 , 
while reproducing the correct value of r T . 

We do not introduce even more delta functions in the low-cr regime, because in such a case we would need values of 
more low-energy QCD observables to fix at least some of the additional parameters. Most of the inclusive low-energy 
QCD observables, with the remarkable exception of the (V + A) ratio r T , either have large experimental uncertainties, 
or large theoretical uncertainties due to large higher-twist contributions, or both. On the other hand, if we fix the 
parameters of the model (with more than two deltas) by simply imposing a further increase in the power index n 
of the difference Ai — a p t ~ (A 2 /Q 2 ) n (n > 6), and without imposing the requirement of the reproduction of any 
additional low-energy observable value, some of the delta's become negative, indicating numerical instabilities. 

In Sec. |TT] we describe the model and impose the conditions which will fix the unknown parameters. In Sec. IIIII 
we explain how to evaluate, in any anQCD model, the leading-twist contribution of the spacelike observables and of 
the related timelike observables, among the latter being r T . In Sec. IIVI we present the numerical determination of 
the model parameters and other numerical results. In Sec. [V] we summarize the results of this work and outline the 
prospects of further applications. 

II. DESCRIPTION OF THE TWO-DELTA MODEL 

As outlined in the Introduction, we construct the two-delta (2d) anQCD model by starting with an ansatz for the 
discontinuity function pi(cr) = Im Ai(Q 2 = — cr — ie) (for a > 0) which agrees with the perturbative counterpart 
pf '(a) = Im a p t(Q 2 = —o — ie) at sufficiently high scales a > Mq (Mq ~ I GeV 2 ); and in the low-scale regime 
< cr < Mq its otherwise unknown behavior is parametrized as a linear combination of (two) delta functions 
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(2d) 



(cr; c 2 ) = nJ2f 2 A 2 6(a - M 2 ) + @(a - M 2 ) x pf\a; c 2 ) (6) 
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where we denoted s = cr/A 2 , Sj = M 2 / A 2 (j — 0,1,2), and rj pt) (s;c 2 ) = pf t \a;c2) = Im a pt (Q 2 = -cr - ie;c 2 ). 
Here, A 2 (~ 10 _1 GeV ) is the Lambert scale appearing in the following expression for a pt : 

I 

~i[l-c 2 /c 2 + W Tl {z)} ' 

where Q 2 = \Q 2 \ exp(i</>), the branches W_i and W+i of the Lambert function refer to the case < co < +n and 



a pt (Q 2 ]c 2 ) = -777; n _iS > w TTZvi ' ( 8 ) 

the bram 
-7r < (f> < 0, respectively, 4 and 
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exp (— iM/c-i) , (9) 



3 In addition to the scale Aqcdi which was fixed by the condition of reproducing the world average value of a p t(M%) = 0.119/7T, in the 
MS scheme. 



In MATHEMATICA ^27j, the functions W n (z) are implemented by the command ProductLog[n, 2] 



where the aforementioned Lambert scale A 2 appears. The explicit expression 
renormalization group equation (RGE) of the form 



is the solution of the (perturbative) 



da p t(Q 2 ;c 2 ) 
<91nQ 2 



-Axipt 



[1 + (ci - (c 2 /ci))a pt ] 



[1 - (c 2 /ci)a 



(10) 
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Here, /3 = (1/4)(11 - 2n//3) and ci = /?i//3 = (1/4)(102 - 38n//3)/(ll - 2n//3) are the universal constants, 
while c 2 = /3 2 //?o is the free three-loop renormalization scheme parameter. The expansion of the beta function 
j3(a p t) = da.pt/ d\nQ 2 in general gives 



/3(a pt ) 



-/3oa pt (l 



cia pt 



C 2flpt 



c 3a P t 



+ ...) 



(11) 



where Cj (j > 2) are general renormalization scheme parameters. The choice of the beta function on the right hand 
side of (|10p gives Cj = c? 2 ~ /cj - (j > 3), which means that the three- loop scheme parameter c 2 can be chosen freely 
in this form, while the higher-loop scheme parameters are then fixed. The specific "effective three-loop" perturbative 
beta function of the Pade form of Eq. (fTU)) was chosen for convenience, because it gives an explicit solution © , in terms 
of the branches of the Lambert function W [28l - t3G |. and, at the same time, it allows us to vary the renormalization 
scheme at the three-loop level (02). In the following, it will turn out to be convenient to vary the scheme parameter 
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c 2 in the explicit solution, the latter being used for p^ 
function @-Q of the anQCD model. 

The Lambert function W — W(z) is defined via the inverse relation z 



Im a pt (Q 2 = — cr — ie) appearing in the discontinuity 
W exp(W), cf. Fig.OJa). The two branches 
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FIG. 1: (a) The defining relation z = We w for the Lambert function W(z), for —l/e < z < 0; (b) The branch W—i(z) for the same 
z-interval; when C2 < 0, the denominator of Eq. (8} becomes zero at a z(«l) in this interval. 



Wipi(z) of the Lambert function are related via complex-conjugation W+i(z*) = W-i(z)*, and the point z = —l/e 
is the branching point of these functions. In the interval —l/e < z < 0, W—i(z) is a decreasing function of z, 
cf. Fig. [Tib). When z -> —0, the scale Q 2 tends to Q 2 -> +cx), cf. Eq. ©, and VF_i(z) — > — oo, this reflecting the 
asymptotic freedom of a p t(Q 2 ) of Eq. ©. In our considered case of low-energy QCD (i.e., with number of quark 
flavors nf — 3), the solution ([SJ has unphysical (Landau) singularities along the positive Q 2 axis, for any c 2 . An 
extension of such beta function to the effective four- and five-loop case, such that the solution is explicit and involving 
Lambert functions, was made in [31(. For more details on the Lambert functions, we refer to Refs. [28l - l32| . 

The aforementioned branching point of nonanalyticity z{si) — —l/e corresponds, according to Eq. @, to the scale 

Q 2 (sl) = A 2 sl with Sl = c\ ° (sl = 0.6347 when nf = 3), and the interval Q 2 € (0, A 2 sl) represents the interval 
of the unphysical (Landau) singularities of a p t (Q 2 ) of Eq. ([8]) . If the scheme parameter c 2 is chosen to be negative (this 
will be our case), then there is an additional pole- type Landau singularity at a somewhat higher scale Q 2 (ui) = A 2 ul 



(<*=> z = z(ul) 



-Po/ci 



/(cie)) at which the denominator of Eq. ([5]) becomes zero, cf. Fig.[TJb), i.e., when 
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•C2/c? = W_i — |« L | 
\c\e 



-/3o/ci 



(12) 



When rif — 3 and c 2 = —4.76 (this will be our central choice of the scheme later), we get ul = 1.0095 (> Sl). For 
this case, the coupling a pt is presented in Fig. [5Ja) as a function of z (for —l/e < z < 0, i.e., slA 2 < Q 2 < oo), and 
in Fig.^b) as a function of t = -ln(-z) = 1.266 ln(Q 2 /A 2 ) + 1.575. 





t=-ln[-z] 



t = 1.266 In + 1.575 

U 2 



FIG. 2: (a) The perturbative coupling a p t of Eq. (O as a function of z, for — 1/e < z < 0; (b) as a function of t = — ln(— z). The curves 
are for the case of n f = 3 and c 2 = -4.76 (=> 4 = 1.266 ln(Q 2 /A 2 ) + 1.575). 

It can be checked that, as a result of application of the Cauchy theorem to the function a p t(Q 2 )/(Q 2 — Q 2 ) in the 
complex-Q 2 plane, the following dispersion relation for a pt holds: 



,rft n 1 f°° , ImOpt(-o-- ie;c 2 ) 1 

a pt (<9 ;c 2 ) = - / ac - 

7T J a =-Ql lin ~ri' 



(a + Q 2 



ds 



s=—u m i n —r] 



r[ pt \s;c 2 ) 
(s + Q 2 /A 2 ) 



(77,7? -+0), (13) 



where the integration covers the entire cut, i.e., starting at a sufficiently low negative value cr min = —Q 2 ^^ (Q^m ~ 
1 GeV 2 ). The perturbative discontinuity function is denoted as r^ p (s;c2) = Im a pt (Q 2 = — sA 2 — ie;oz). Since the 
cut of the coupling a pt (Q 2 , c 2 ), Eq. ([8]) with C2 < 0, includes also the pole Q 2 = ulA 2 of the coupling, the contour 
of integration in the complex (Q 2 /A 2 )-plane is of the type as presented in Fig. [3] (with the outer radius going to 
infinity). Therefore, the dispersive relation (J13I) obtains a slightly generalized form 
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FIG. 3: The path of contour integration in (Q 2 /A 2 )-plane leading to the expression (TH)) . 



a P t(Q 2 ) = - 



*^S + ^r *- 



a pt ((u L + Aue^)A 2 ;c 2 ) 

i$l 



ds 



(s + Q 2 /A 2 ) • 2tt y $= _, " ' ■ [(g 2 /A 2 ) - u L - Aue 

r { ^\s;c 2 ) Res {z=UL) a p t(zA 2 ;c 2 ) 

{s + Q 2 /A 2 ) + (- UL + Q2/ A 2) • 



(14) 
(15) 



In Eq. (|14p . AuA 2 is a sufficiently small (but otherwise arbitrary) radius of integration around the point ulA 2 in the 
complex Q 2 -plane, cf. Fig. [3J In Eq. (fTS")) . this integration is expressed by the residue of the function a pt (zA 2 ) at this 
point. 

The perturbative discontinuity function r\ p '(s; c 2 ) = Im a p t(Q 2 = — sA 2 — ie; c 2 ), which is nonzero for — sl < s < 



-oo, has the specific form 
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The analytic (spacelike) coupling Ai(Q 2 ;c 2 ) of the two-delta anQCD model is constructed on the basis of the 
discontinuity function (J7J [cf. Eq. (ITBl ] using the dispersion relation (JT|). This gives 



A 1 (Q 2 ;c 2 ) = J2 



fj 



i (U + Sj) 



„(pt) 



(is 



(s;c 2 ) 



(s + u) 



(18) 



where w = Q 2 / A 2 . 

In the presented two-delta anQCD model, we will consider the first three quark flavors to be massless, and will 
consider that the momenta of the rij — 3 regime in the anQCD model reach up to the threshold \Q 2 \ = (2m c ) 2 
(ps 6.45 GeV 2 ). Further, the anQCD model will be constructed in such a way as to practically merge with pQCD at 
such sufficiently high momenta [in the considered renormalization scheme as fixed by Eq. (|10[) ]. Therefore, we will 
consider that the value of the Lambert scale A 2 used in our analytic coupling A\ (Q 2 \ c 2 ) coincides with the perturbative 



Lambert scale A pt , the latter being determined by the condition 



4" s) (^i 



0.1184/7T, i.e., by the central value 

(MS), 



of the world average [33j. Therefore, A 2 is determined by RGE-evolving this a pt down to a^ ((2m c ) ;n/ = 3), 

using the four- loop polynomial form of /3' MS )(a p t), and the three-loop matching conditions [34| at quark thresholds 
/i 2 = (2m q ) 2 (q = b, c); and then changing from the MS scheme to the scheme 'c 2 ' [= (c 2 , c^ — c 2 /ci, . . .)} defined by 
the beta function of Eq. (fT0|) (as explained, e.g., in Refs. [18l.[26l|V 
The conditions we impose to fix the parameters are the following: 

1. The analytic coupling Ai(Q 2 ; c 2 ) acquires the aforementioned pQCD value of the scale A 2 of a p t(Q 2 ; c 2 ;n/ = 3) 



at Q 2 = {2m c y 



A 2 =A 2 pt (n f = 3) 



(19) 



2. While in general we expect Ai(Q ;c 2 ) to differ from a pt (Q ;c 2 ) at Q > A 2 by ~ (A 2 /Q ) , we impose the 
condition 



A(0 2 ;c 2 )-a pt (Q 2 ;c 2 )^(A 2 /Q 2 )' 



with 



= 5 . 



(20) 



The condition (f2"0")) represents in practice four conditions. In the following few lines we will describe how to formulate 
these four conditions. 

When we subtract from the perturbative coupling (fl5|) the analytic coupling (fP8)l , we obtain 



a pt (Q 2 ;c 2 )-.4i(Q 2 ;c 2 ) 



Res 



(z=ul) 



a pt (zA 2 ) 



(u - ul) 
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/ 



\ {U + Sj) 



ds 



r\ (s;c 2 ) 
(s + u) 



(21) 



Expanding the left-hand side in powers of (l/u) = (A 2 /Q 2 ), the imposition of the condition ([20]) gives us the conditions 
that the terms of ~ (A 2 /Q 2 ) 1+fe (fc = 0, 1, 2, 3,) in this expansion give us zero, i.e., we have the following four conditions 






d SS M pt) ( S ;c 2 ) + (- UL ) fe Res (z=UL) a pt (zA 2 ) = «i/i + 4fl (k = 0,1,2,3) 



(22) 



Altogether, Eqs. (fTi?]) and ((2^|) represent five conditions. Once the scheme c 2 parameter is chosen, we have altogether 
six parameters in the model: / 2 , /f, si, s 2 , sq, and the scale A. Therefore, yet another condition will have to be 
imposed, possibly involving the correct reproduction of a low-energy observable. The best candidate for this appears 



to be the canonical (V + A) nonstrange and massless semihadronic r-lepton decay ratio r T , (35I437J ]. 5 When we remove 
the (measured) strangeness-changing contribution, the color and CKM factors and the electroweak effects, and the 
chirality- violating higher twist (quark mass) contributions, the following value is obtained (cf. |38|-|4(j], and Appendix 
B of Ref. |l| for details) 

r T (AS = 0,m ? = 0) cxp . = 0.203 ± 0.004 . (23) 

Numerical analyses of the measured data indicate that the chirality-conserving higher-twist effects, such as gluon- 
condensate contributions, are negligible in the case of the considered V + A decay channel. Although such analyses 
have been performed within pQCD+OPE approach, we will assume that they remain valid when an analysis is 
performed with the presented anQCD two-delta model + OPE. This assumption appears to be reasonable because 
the considered anQCD coupling Ai(Q 2 ) is very close to the pQCD coupling a p t(Q 2 ) (in the considered scheme) at 

momenta \Q 2 \ ~ 1 GeV 2 . Therefore, in the calculation of the discussed r T , Eq. (|23|) . in the presented anQCD model, 
the (chirality-conserving) higher-twist contributions will be ignored. 

III. CALCULATION OF ADLER FUNCTION AND r T IN ANALYTIC QCD 

The calculation of r T is then performed in the same way as presented in |18l. |26J|. i.e., by performing explicitly the 
integration corresponding to the leading- /3 (LB) resummation for r T , and adding the three known beyond-LB (bLB) 
terms (i.e., including the exact contributions of ~ ■A&). 

In this Section we will present only the main points of calculation of spacelike (such as Adler function) and timelike 
quantities (such as r T ) in anQCD models. For details, we refer to [ill, [lj, [26| , and especially Appendices C and D of 

m 

The basic idea of the approach in the evaluation of spacelike observables T>(Q 2 ) in general anQCD model is to 
reorganize first the perturbation series r D(Q 2 )- pt 

£>(<9 2 ) P t = a P t + diapt + a^apt + • • • , (24) 

into the modified perturbation series (mpt) 

£>(Q 2 )m P t = a pt + dia p t,2 + c?2a P t,3 + ■ ■ • , (25) 

where dfca P t,fc+i are the logarithmic derivatives of a p t 



^ (-l) k d k a pt (Q 2 ) 
$k\ <9(lnQ 2 )* 



S P t, fc+ x(Q 2 ) ^ ^ U J^> , (fc-0,1,2,...) , (26) 



They are related with the powers of a p t = a s /n in the following way [using renormalization group equations (RGE) 
inpQCD]: 

Spt,2 = «pt + cia-lt + c 2flpt + ■ • ■ » (27) 

5 
a pt , 3 = ap t + -Cia pt + . . . , a pM = a* t + . . . , etc. (28) 

This, in turn, means that the powers of a p t are linear combinations of logarithmic derivatives 

a pt = a pt , 2 - CiSpt,3 + ( ^c 2 - c 2 J a pt , 4 + ■ • ■ , (29) 

5 
Opt = a pt , 3 - 2 C i fl pt,4 + ■ • • , a^ = a pM + . . . , etc. , (30) 

which allows us to relate the mpt coefficients with the usual perturbation series (pt) coefficients 

di = di , di = di — c\d\ , (31) 

d 3 = d 3 - -c x d 2 + f -c 2 - c 2 J d\ , etc. (32) 



5 The "canonical" means that the normalization is used such that (r T ) p t = a p t + 0(a? >t ). 
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In [ll|, [I4I it was shown that it is imperative to construct first the analogs of the logarithmic derivatives of a pt in the 
following way: 6 

I 0(lnQ')* J an = 9(lnW ( Pt(0 ))a " = 0(hW (fc = 0,l,2,...). (33) 

This means that the 'mpt' expansion (P25l) becomes in anQCD the corresponding "modified analytic series" ('man') 

2?(Q 2 )man = A x + d x A 2 + d 2 A 3 + . . . , (34) 

and its truncated version is 

V(Q 2 )lnL = Ai + diA 2 + ... + d N -!A N . (35) 

Here we denoted by Ak+i the logarithmic derivatives of A\ 

~ 2 (-l) fc a%(M 2 ) ,, 1 . , ,_ 

A +1 (m) = ^ d{lnfi2)k , (fc = l,2,...). (36) 

The expressions (l54"l) - ([3"o| are the basis of the evaluation of massless spacelike observables in any anQCD. 7 

The quantity r T is, on the other hand, a timelike observable. However, it is expressed via a contour integration [41( 

r r = ^J n d<P(l + e l0 ) 3 (l - e*+) d Adl (Q 2 = m 2 .^) , (37) 

which involves the (spacelike, massless, nj = 3) Adler function dAd\(Q 2 ) — Opt(Q 2 ) + C( a P t)- The perturbation 
expansion of dAdi is known up to ~ <Zp t |42h44| | 

3 

dAdi(Q 2 ) p 1 = Opt + ^(dA dl )„a;\ +1 • (38) 

71=1 

On the other hand, the leading-/3 parts (dAdi)n = Ch,nPo [= (^Adi)n ] of all the coefficients (dAdi)n are known 8 
[45l l46l ] and the resummation of these contributions can be performed by using the approach of Neubert of integration 
with characteristic functions [47| - this can be performed in any anQCD without ambiguities (since no Landau 
singularities) 

(dAdiffi (Q 2 ) = A, (Q 2 ) + c[%MQ 2 ) + ... + c^A n+1 (Q 2 ) + ... (39) 

= / -F A di(tMi(i0 2 e- 5/3 ) . (40) 

Jo * 

Here, i*Adi(£) is the characteristic function of the Adler function, whose explicit expression was obtained in |47| . 

On the other hand, the full coefficients (e?Adi)n and (dAdi)n are known only up to n = 3 [42T - l44j |. Therefore, the full 
Adler function can be evaluated in anQCD by adding to the leading-/? (LB) contribution the three known beyond-LB 
(bLB) terms 

(d A di)£^ bLB) (Q 2 ) [4] - / -f\, u li)A^i()', ' ! ) > ;iA (1 il„.4„-, . Nil 



/ HE Ff dl (t) A(tQ 2 e" 5 / 3 ) + V(T Ad i)„^+i , 

J J n=l 



6 If the analytization is performed in any other way, the renormalization scale and scheme dependence of the resulting truncated analytic 
series of any observable T>(Q 2 ) will in general increase (instead of decrease) when the number of terms in the series increases, cf. [Ill,ll2fl. 

7 If masses are involved in the evaluated physical quantity, perturbation series contains noninteger powers a" t (and possibly derivatives 
thereof with respect to v, i.e., a" t ln k a p t)- The evaluation of such quantities in anQCD models is then based on the procedure presented 
infJJ. 

8 We have (rfAdl)n = c„ : „/3q + 0(/3q~ ) and (d\di)n = (ii,nPo + ^(Pq~ )• The expansions in powers of /3o are obtained when (d^dOn 
and (c(A.dl)n are written in powers of nj (= — 6/3q + 33/2) and then reorganized in powers of /3q. 



where 

(T A dl)„ = (dAdl)n ~ c$0» (42) 

are the complete bLB parts (~ /3g _1 ) of the coefhcients (dAdi)n- 

Using the expression (1411) in the contour integration (J37I) gives for r T 

r ^B+bLB) lW = r (LB) + £ (TAdl ) n7( l n+l) , (43) 

ra=l 

where 7(„4 ra +i) are the contour integrals given by 

I(A n+ i) = ^- I *c#(l + e*) 3 (l-e**);Wi(m?e*) (n = 1,2,3), (44) 

and the LB part in (J43J) is a well-defined (in anQCD models) integral of the form 

r (LB) = I r^F r (t) Pl (tm 2 T e- 5 / 3 ), (45) 

where F r (t) was calculated in [18| from the Minkowskian characteristic function F^ 4 (t) of [48| . For details on F r (t) 
we refer to App. D of [18(. 9 

Thus, the more explicit form of the r r -reproduction condition, mentioned in the previous section, is 

^(LB+bLBW) = ) r (LB) + j^ (T Adl ) n I(A n+1 ) = 0.203 . (46) 

n=i 

The six conditions, (IT9"|) . (j2"2")l and (BUI) , then determine the six parameters of the model: /?, /|, Si, S2, so and the scale 
A. This procedure can be performed once we have chosen a value of the scheme parameter c^ of Eq. (ITD")) . 

IV. NUMERICAL RESULTS 

In this Section we present the numerical results for the parameters of the model, obtained from the imposition of 
the aforementioned six conditions. The additional (implicit) conditions that we choose are that the weights ff and 
/f are positive. This is based on the fact that the discontinuity function pi(a) is positive in any reasonable scheme of 
pQCD. Furthermore, the condition of positivity of pi(cr) can be expressed also via the condition that the Minkowskian 
coupling 

8 lWB if^ (47) 

7T Jo a 

is a monotonously decreasing function of scale a. 

If we choose for the scheme parameter C2 = 0, we obtain from the six conditions, as a result, that the pQCD-onset 
mass Mq = -ysiJA is relatively high, Mq ~ 1.7 GeV. Increasing C2 to positive values increases Mq further; e.g., for 

c 2 = c^ s /2 (sa 2.24) we get Mo ~ 2.0 GeV, etc. On the other hand, decreasing c 2 to negative values, we obtain 
smaller Mo. We believe that the effective pQCD-onset scale Ma should be significantly smaller than the mass of the 
r lepton. This turned out to be so in the one-delta anQCD model of (26(, with c 2 = 0, where Mo &w 1 GeV was 
obtained. In general, we do not want to parametrize (via delta functions) relatively small deviations from pQCD, i.e., 
those at \fa > 1.5 GeV. On the other hand, by lowering the value of C2, we encounter at C2 < —8.0 (Mo < 0.9 GeV) 
negative values of S2, implying that the analyticity is lost. 

Therefore, we will adjust the scheme parameter C2 in such a way as to get Mo = 1.25 ± 0.25 GeV. The results are 
given in Table |U We believe that all three choices of M will give almost the same predictions for various physical 



9 In I|43j l- ll44ll we used in the Adler function the renormalization scale fi 2 = Q 2 (= m^e 1 * 1 ); this scale, of course, can be varied, cf. [18 
for details. 
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TABLE I: The parameters of the considered two-delta anQCD model, for the three chosen values of the pQCD-onset scale Mo: 
1.0; 1.25; 1.50 GeV. In addition, the results of the one-delta (Id) model of [2g] in the scheme c-i = are given (the last line). 
See the text for details. 



M [GeV] 


C2 = 


<h/f3o 


A [GeV] 


so 


Si 


f? 


S'2 


/! 


Ai(0) 


1.00 
1.25 
1.50 


-7.15 
-4.76 
-2.10 




0.193 
0.260 
0.363 


26.86 
23.06 
17.09 


19.473 
16.837 
12.523 


0.3637 
0.2713 
0.1815 


0.3594 
0.8077 
0.7796 


0.7808 
0.5409 
0.3462 


2.29 

0.776 

0.544 


0.886 (Id) 







0.472 


3.525 


0.4755 


0.2086 






0.544 



observables. The reason for this lies in the renormalization scheme independence of pQCD results; and our model, 
although nonperturbative and analytic at low momenta \Q\ ~ 1 GeV, is practically indistinguishable from pQCD 
at all higher momenta. It is interesting that the value of the coupling A\ at Q 2 = 0, obtained in this model for 
M = 1.25 GeV (-4i(0) ~ 0.8) is not far from the value in Ref. [3] (Ai{0) ~ 8.9/N c /n sa 0.9), which was obtained 
from an analysis using Dyson-Schwinger equations for the ghost and gluon sector under an assumption of regularity 
of the ghost-gluon vertex. 

In Table U we included, in the last line, also the results of the one-delta (Id) model of [26|], obtained in an analogous 
way, in the scheme c% = C3 = . . . = 0. 10 In that model, though, the smaller number of parameters led to less stringent 
conditions (J2T)]) . namaly with n max = 3 (cf. also a similar model in [50]). 

In Fig.HJa), we present, for the resulting "central" choice of M = 1.25 GeV (c 2 = —4.76), the corresponding pQCD 
discontinuity function pf '(a) = Im a pt (— a — it) of the underlying perturbative coupling Opt ©, in the regime of low 
|cr| and including the unphysical (Landau) regime of negative-cr cut. We can note that in the latter regime there is an 
additional, pole-like singularity at a = — i*lA 2 (« —0.261 GeV ); while the "continuous" part of pf (c) ends a bit 
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FIG. 4: The discontinuity function pi(cr): (a) for the underlying pQCD coupling a p t(Q 2 ), where the Landau pole at a = — ulA 2 
(as —0.068 GeV 2 ) and the branching point at a = — slA 2 (fs —0.043 GeV 2 ) are visible; (b) for the considered two-delta anQCD model, 
where the unknown region < a < Mq is parametrized by two delta functions, at a = M 2 and a = M^- The parameters used correspond 
to the central case in Table[I](Mo = 1-25 GeV). 



earlier, at a = — s^A 2 w —0.207 GeV 2 . 11 In Fig. SJb), pi(cr) of the considered two-delta anQCD model is presented, 
cfEq. @. 

In Figs. [Ha) and (b), we present, for the aforementioned central case Mq = 1.25 GeV, the resulting spacelike 
coupling Ai(Q 2 ) at positive Q 2 . The higher order couplings Ak{Q 2 ) (k = 2,3), cf. Eq. (|36)l . are also presented; due 
to a strong hierarchy (even at low Q 2 ) they are rescaled, for better visibility, by factors 4 and 16, respectively. All 
the corresponding pQCD quantities (a p t, 4a p t,2, 16a p t,3) are presented as dotted curves. For better visibility at low 
Q 2 , Fig.(SJb) is presented with Q 2 on logarithmic scale. It is clearly visible that the model practically agrees with the 



10 In Ref. Lffl . the parameters differ a little from those in Table |I] because the world average value taken there was from the year 
2008, Opt(Mi; MS) = 0.1190/7T, [49| : and because we imposed there the condition _4i((3m c ) ) = o p t((3m c ) \ n f = 3) instead of the 
(numerically similar) condition (| 196 . 

11 Note: sl = C\ C1 ' , which is approximately 0.635 when rif 
this is the branching point for a p t(Q 2 ), cf. [281. l3ll] ■ 



3; at Q 2 = slA 2 the Lambert function W-i(z(Q 2 )) is equal to —1, and 
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FIG. 5: The analytic couplings Ai (full line), 4 X A2 (dashed curve), and 16 X A3 (dash-dotted curve) for positive Q 2 , in the considered 
two-delta model, for the central parameter choice of Table U (Mo = 1.25 GeV): (a) linear scale is used for Q 2 ; (b) logarithmic scale is 
used for Q 2 , for better visibility at low values of Q 2 . For comparison, the corresponding pQCD couplings (a p t, 4 X a p t,2, 16 X a p t,3) 
are included, as dotted curves (in the same renormalization scheme, with C2 = —4.76); the vertical dotted line is the Landau pole at 
Q 2 = Ul A 2 (« 0.068 GeV 2 ). 

corresponding pQCD model at Q 2 > 1 GeV ; and that for Q 2 < 1 GeV the model differs from pQCD significantly, 
due to the imposition of the analyticity. 

In Fig. EJa), we present the difference between the perturbative and the analytic coupling (a pt (Q 2 ) — A\(Q 2 )) at 
positive Q 2 , for the central choice of parameter (Mo = 1.25 GeV) of Tabic Q] In the Figure we keep, formally, n/ = 3 
even for high Q 2 . We see that the difference vanishes fast when Q 2 increases. For comparison, in Fig. Efb) we show 
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FIG. 6: The difference between the perturbative and the analytic coupling, as a function of positive Q 2 , (a) the considered two-delta 
anQCD model (with Mq = 1.25 GeV); (b) for the one-delta anQCD parameter of Ref. ;26i] , but with the input parameters somewhat 
modified (see the text for details). 



the difference (o p t — A\) in the case of the one-delta anQCD model of 26], with parameters as given in the last line 
of Tabled] 

From Fig.[Sfa), we can deduce that the difference (a pt — A\) behaves as oc (A 2 /Q 2 )™ oft , where the numerical value 
of the effective power index n e g is somewhat less than 5; it is somewhere between 4 and 5. This can be understood 
in the following way: the values of sq and si in the model are relatively large (~ 10 1 ), and therefore the coefficients 
K, n (for n — 5, 6, ...) in the expansion 



a v t{Q ; 02) ~ MQ 2 ; c 2 ) = IC 5 (A 2 /Q 2 ) b + /C 6 (A7QT + 



(48) 



are large (~ s$~ , s™ _ ) and are increasing when n increases. 

Finally, in Table [TT| we compare the convergence of the series for r T , Eq. (J4U]) , in the three cases of the two-delta 
model and in the one-delta model of Table fl] We see that the convergence is good in general, and appears to be better 
when the underlying pQCD scheme parameter C2 is more negative. 

There is one more aspect of the presented model that we should address. We recall in more detail how we matched 
the presented anQCD model indirectly to the MS scheme, via the condition (JTI 
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TABLE II: The four terms in the truncated analytic expansion (|43[) in two-delta anQCD model for the three choices of the 
pQCD-onset scale Mo and for the one-delta anQCD model, as specified in Table H] 



Mo [GeV] 


C2 


A [GeV] 


r T : LB 


NLB 


N 2 LB 


N 3 LB 


sum (sum) 


1.00 
1.25 
1.50 


-7.15 
-4.76 
-2.10 


0.193 
0.260 
0.363 


0.2023 
0.2052 
0.2096 


0.0012 
0.0013 
0.0013 


-0.0010 
-0.0053 
-0.0110 


0.0005 
0.0019 
0.0030 


0.2030 
0.2030 
0.2030 


0.886 





0.472 


0.2149 


0.0014 


-0.0176 


0.0043 


0.2030 



1. The Lambert scale of the model, at n/ = 3, was chosen to coincide with the Lambert scale of the underlying 
pQCD model in the same scheme determined by the beta function of Eq. (fTUf at nf = 3; 

2. The latter scale was fixed in such a way that the change of the scheme parameters (c 2 , c 2 /ci, c 2 /cf , . . .) to the 



four- loop MS scheme parameters (02,03), at n/ = 3 and fi 



(2m c ) 2 , gave the value of Op t ' ((2m c 



n f = 3) 

; 0.1184/-7T, the central value of the world average [33J]. The RGE-running 
between /i 2 — M'% and /j 2 = (2m c ) 2 , in MS, was performed in the usual way, using the four-loop polynomial 
form of /3 < - MS ' ) (apt), and the three-loop matching conditions [34| at quark thresholds /1 2 = (2m q ) 2 (q ~ b,c). 



which corresponds to a pt (M§) 
Ml and \i 2 



As pointed out in Ref. [511 ] . the aforementioned (three- loop) matching in principle introduces, indirectly, an element 
of nonanalyticity in the described framework, at 12 the scale Q 2 — (2m c ) 2 and (to a much lesser degree) at Q 2 = (2ml,) 2 . 

This is so because the matching introduces nonanalyticity (even: discontinuity) in the running coupling a t (Q 2 ) 
at those threshold scales. It would be more convenient for the presented low-energy anQCD model to be matched 
[at energies Q 2 ~ (2to c ) 2 ] to a scheme which introduces the quark mass threshold effects in the running in a gradual 
(analytic) way. One such scheme is the pinch technique (PT), Refs. I52| . The PT effective charge (i.e., running 
coupling) was presented for supersymmetric QCD at one-loop in Ref. [5l| . Specifically, in nonsupersymmetric QCD 
with three massless flavors, the relations of Ref. [5l| (especially their Eq. (A3)) imply (for |Q§|, |Q 2 | < (2m c ) 2 ) 



a PT (Q 2 ) = a PT (Q 2 ) 



t PT 



(Qo) 



\hQ 2 /QD + \li/2{Q 2 I 



-L 1/2 (Q 2 /m 2 c ) 



where the quark threshold function L 1 i 2 is 

L 1/2 (Q 2 /m 2 ) = (3 - fi 2 ) [fi ArcTanh(l/^) - 1] + 2 



(49) 



(50) 



where /3 = \J\ + Am 2 /Q 2 and the (complex) momenta Q 2 are not on the cut Q 2 < —Am 2 . In the case of decoupling 
(|<2 2 | <C rn 2 ), this function acquires the value of (5/3). 

We can estimate the practical errors introduced in the calculations in our model due to the aforementioned matching 
to MS, by estimating the errors that the actual nondecoupling of the charm quark mass introduced in our calculation 
of r T (note that in our calculation, we considered c quark to be completely decoupled, in the spirit of MS). The 
relation (|49j) would be transcribed in our model, approximately, in the following way: 



-4i(Q 2 )thr « -4i(Qg)thr + A 2 {Q 2 ) 



1 



1 



HQVQo) + 7 L 1/2 {Q 2 /m 2 c ) - -L 1/2 (Q 2 /m 2 c ) 



(51) 



where A 2 (Q 2 ) is the logarithmic derivative of A\ defined via Eq. (|36]h and the scale Q 2 is taken to be the geometric 
mean 13 of Q\ and Q 2 : Q 2 = \/Q 2 Q 2 . Further, the subscript "thr" indicates that the nondecoupling effect m c 7^ 00 is 
taken into account, in a first approximation. The corresponding relation in the model without the analytic threshold 
effects is 



-4i(Q 2 ) « AxlQl) + A 2 (Q 2 ) 



-\hq 2 /q 2 ) 



(52) 



12 We consider, throughout, the first three flavors to be massless. 

13 Formally, at one-loop level, it would be equivalent to take the scale Q 2 = Qq or Q 2 = 
these two scales (i.e., the arithmetic mean of the logarithms of these scales) is better. 



Q ; numerically, though, the geometric mean of 
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TABLE III: The estimates of nondecoupling of m c in the calculation of r T , Eq. (|54|) . for different value of the parameter n. In 
the third column we include a measure of the effect of nondecoupling of m c at Q\ — smj (assumed to be negligible in Eq. (|54l) ); 
and in the fourth column the ratio of | ln(Q 2 /Qo)A2{Q 2 ) with Ai(Qo) (assumed to be small in (|51[1 and (|52p ). for Q 2 — m 2 
and Q — «%■ 



K 


Sr^°\m c7 t 


oo) 


(L 1/2 (Kml/m 2 c ) - 


-L 1/2 (0))/L 1/2 (0) 


- 1 ln(K)^l 2 {-jKml)/A\ (Km 2 .) 


0.15 


6.1 x 10 -4 




0.034 




0.42 


0.20 


4.6 x 10 -4 




0.045 




0.36 


0.30 


3.1 x 10" 4 




0.066 




0.28 


0.40 


2.3 x 10 -4 




0.087 




0.21 



An estimate of the error in the calculation of r r , due to the nondecoupling of m c , in the calculation of r T can be 
obtained by using the relation (|37]) in the leading order (LO) when dAd\(Q 2 ) = A\(Q 2 ) there, and subtracting for 
A\(Q 2 ) the two relations ([51]) and (j52")l on the contour 



5r^°\m c ^oo) = i- J ^0(l + e^) 3 (l-e^) [A(Q 2 )thr - -4i(Q 2 )] 

d0(l + e^) 3 (l-e^) 



1 
2^ 



(A(Qo)thx - -Ai(Qo)) + g-42(VQ^Q 2 ) (L 1/2 (Q 2 /m 2 c ) - L 1/2 (Q 2 /m 2 c )) 



\Q 2 =m 2 e 



.*-(53) 



The scale Qo nere should have a low value (|Qg| <C m 2 ) so that the charm quark can be assumed to basically decouple 
(m c = oo) at such scale. Namely, in such a case, the first (unknown) term in the integrand, (-Ai(Qo)thr — Ai(Qq)), 
can be taken to be zero, and the remaining term would indicate correctly the effects of nondecoupling of m c . On 
the other hand, |Qo| 2 cannot be taken to be too small, as then the ~ A 2 term in the relations (I5T1) and ([5"Hj) would 
dominate over the A\{Q%) term, due to a very large value of the logarithm \n(Q 2 /Q 2 ). Therefore, we will choose 
Qq = kQ 2 (= Km 2 , exp(z^)), with k < 1 varying in a specific interval were both aforementioned restrictions are 
reasonably fulfilled 



Sr^°\m c 



>"»»"i5F 



+ 7T 



(1 + e^) 3 (l - e^)A 2 {jQ 2 Q 2 ) (L 1/2 (Q 2 /m 2 c ) - L 1/2 (Q 2 /m 2 c )) 



Q 2 



b ;Qh 



(54) 
The results of these estimates, for the central value of the parameters of the presented two-delta anQCD model 
(i.e., the case Mo = 1.25 GeV of Table H]) are presented in Tabic UTTl We see from the Table that the effects of the 
introduction of analytic threshold effects (i.e., the effects of the nondecoupling of m c ) in the calculation of r T are 

Sr T < 10~ 3 . This is to be compared with the theoretical value r r w 0.123 in the anQCD model, and the full 
value r T — 0.203 in the model, and with the experimental uncertainties 8r T — ±0.004 [Eq. ([23]) ]. We conclude that, 
although these effects are appreciable, they are somewhat lower than the present experimental uncertainty of r T . 



SUMMARY 



In this work we presented an analytic QCD (anQCD) model which, at high squared momenta \Q 2 \, becomes to 
a high degree indistinguishable from perturbative QCD (pQCD); nonetheless, at low \Q 2 \ the spacelike couplings 
A. n {Q 2 ) of the model mirror correctly the analytic properties of the (to-be-evaluated) spacelike observables T>(Q 2 ) as 
dictated by the general principles of the local quantum field theories. The model reproduces correctly the experimental 
value of the r lepton (nonstrange, V + A) semihadronic decay ratio r T , i.e., the only low-momentum QCD observable 
that is well measured and whose higher-twist effects arc small and appear to be under control. The difference between 
the analytic Ai(Q 2 ) and its perturbative counterpart a pt (Q 2 ) [= a s (Q 2 )/Tr] is formally ~ (A 2 /Q 2 ) 5 at \Q 2 \ > A 2 , 
where A 2 ~ 0.1 GeV is the QCD (or: light meson) scale. The starting point was the construction of the discontinuity 
function pi(cr) for the analytic coupling A\(Q 2 ). This pi(a) is assumed to coincide with its pQCD counterpart for a 
above a pQCD-onset scale Mq ~ 1 GeV 2 . On the other hand, the unknown behavior of pi(<r) in the low-c regime is 
parametrized as a linear combination of two delta functions. The underlying pQCD scheme parameter c 2 is adjusted 
so that the pQCD-onset scale M is either 1.0, 1.25, or 1.50 GeV. We believe that these three variants of the two-delta 
anQCD model represent almost the same physics, since they reproduce the same value of the key low-momentum 
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observable r T (V + A) and, at the same time, they differ in the value of the underlying pQCD (i.e., high- momentum) 
renormalization scheme parameter C2, the latter being the only free input parameter of the model. If the value of C2 is 
further adjusted so that Mo is below 1.0 GeV, the analyticity gets lost at Mq « 0.9 GeV. The model is an extension 
of the one-delta anQCD model of [26| where c-i — was used. 

The main motivation behind the construction of the considered anQCD model is that with it we can eventually 
evaluate low-momentum QCD quantities whose higher-twist contributions are appreciable. An example are the 
separate vector (V) and axial vector (A) channel of r T . Such evaluations would involve the considered anQCD model 
together with the Operator Product Expansion (OPE). This application (anQCD+OPE) is consistent, due to the very 
suppressed difference Ai(Q 2 ) — a pt {Q 2 ) ~ (A 2 /Q 2 ) 5 in the ultraviolet regime; the latter implying that the higher-twist 
terms in OPE, of dimension d < 8, will still be of infrared origin, in accordance with the philosophy of the ITEP 
(Institute of Theoretical and Experimental Physics) group [23[ . 
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